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Abstract 

The evolution this article is dramatic. In order to estimate the average speed 
of mosquitoes, a simple experiment was designed by Richard (Lu-Hsing Tsai), 
Tom (Po-Yu Tsai) and Robert (Hung-Ming Tsai). The result of the experi- 
ment was posted in the science exhibitions Taichung Taiwan 1993. The average 
speed of mosquitoes is inferred by the simple relation m — KpcAt , where m 
is the number of mosquitoes that are caught by the trap, p is the density of 
the mosquitoes in the experimental device ( a camp tent), c is the mean speed 
of the mosquitoes, A is the effect area of the trap, t is the time interval, and 
A" is a proportional constant. In order to find the value of K. The identity 
E(l) = aV/A is discovered, where V is the volume of a bounded region in the 
space R n , A is the boundary area of the region, a is a proportional constant, 
E(l) is the expectation of the random paths. By a random path in the region 
we mean the path which is traced by a random walker who passes through the 
region. In this paper, we will show how to get the data of a generated by com- 
puter. Actually, a = 1/K. Though the rigorous proof is not shown, a sketch 
proof will be shown in this paper. The theoretical values of K are 1/2, 1/tt, 
1/4, 2/(37r), — , for dimension n = 1, 2, 3, 4, 

Introduction 

Although random walk is just a simple mathematical model, it can be used 
to describe many physical phenomena such as Brownian motion and diffusion. 
There is a very important number called "Boltzmann constant" in physics. His- 
torically, scientists used the random walk model to evaluate this constant for 
the first time [4] . The Bertrand paradox can be traced back to the Buffon nee- 
dle problem established in 1777 [1]. Since then, the research of Buffon needle 
problem has been developed into a branch of mathematics, called geometric 
probability. The most important subject in geometric probability is "random 
chords" [1]. It is very difficult to design an experiment based on its original 
definition because an infinite long straight line is required. To the best of our 
knowledge, no one has mentioned the relationship between the random walk 
model and the random chord problem. In this article, we use a simple way to 
show that random walk problem is essentially a random chords problem. 

A simple experiment 

In 1993, Hung-Ming Tsai [6] had an idea to estimate the mean speed of 
mosquitoes. The idea was: mosquitoes were put in a hemispherical tent with 
a diameter of 230 cm, and a trap was set up somewhere in the tent. When 
a mosquito touched the trap, it would get killed because there we electrified 
grids on the trap. And a electric spark made a sound which was sensed by 
human ears. After the mosquitoes were put in the tent uniformly, the switch of 
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the trap was turned on, the number of sparks was counted, and the time was 
recorded. The idea was nothing other than an experiment. By this experiment, 
he estimated the mean speed of mosquitoes. 

What described above is a crude but interesting experiment. Although the 
experiment is easy, the mathematics for analyzing the data obtained from the 
experiment seems not trivial. Let t be the time interval in which there are m 
mosquitoes that are trapped. Let p be the density of mosquitoes in the tent, 
let c be the mean speed of the mosquitoes and let A be the effective area of the 
trap. In order to analyze the experimental data, we assume 



where K is a proportional constant. From the experiment, we can easily get 
data for m, p, A, and t. If we assume that the flying of mosquitoes is a three 
dimensional random walk, by applying the theory of geometric probability to 
compute the value of K , then the mean speed of mosquitoes can be estimated 
from equation (fT]). We find that it is nontrivial to find the theoretical values of 
K . Due to the Bertrand's paradox and the difficulty, we compute the value of 
K by computer simulation 

The value of K and a new mathematical identity 

Let E(l)be the expectation of the length of random paths in the bounded 
region. Without loss the generality, we assume the bounded region is a sphere. 
In order to simplify the sketch proof of the identity, we use discrete model. If it 
is necessary to treat the continuous case, then the radius of sphere is taken by 
limit process. Let r = E(l)/c , where c is the mean speed of mosquitoes. Then 
t is the mean time interval that a mosquito stays in the sphere. Consider the 
sequence of time to , t\ , . . . , U, . . . , where ti = i ■ r , i = 0, 1, 2, Imagine 
that at each time U , there are exactly N (new) mosquitoes get into the sphere 
and N (old) mosquitoes get out of the sphere. So, at any time, there are exactly 
N mosquitoes in the sphere. By the definition of density, we have p — N/V 
, where V is the volume of the sphere. Now, we consider that the surface of 
the sphere is the sensor whose area is A. During any time interval of length r , 
there are exactly N mosquitoes passing through the boundary and getting into 
the sphere. So the boundary ( or the sensor ) senses Amosquitoes, during a 
time interval of r, that get into the sphere and hence senses Nt/r mosquitoes 
in any time interval of t. Since Nt/ris the m in equation (jTj) , p = N/V, and 
t = E(l)/c, equation (fT]) becomes 



The mathematical identity is obtained so easily that we can not almost 
believe it. Now, we can prove the identity. Since the system 
is in equilibrium state, the number of mosquitoes in the bounded 
region is not changed in any time. Therefore, Am — AN . 

By so called dimension analysis, we have three freedoms to choose 
the scale in equation (fT]) (a) We can take the time scale so that 



m = KpcAt 



(1) 



K = A/{VE{1)) 



(2) 
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there are N mosquitoes walk into and walk out the bounded region 

in a unit time, that is, N — KpcA • 1 . (b) We can choose the speed of random 

walker is 1, that is, c = 1 in random walk, (c) By the definition of density, 

p = n/V the quantity E(l)A/V is dimensionless. We can choose 

the unit of the length so that E{1) = 1 and hence E(l) = c. Therefore, 

equation (J) becomes 1 = KAE(l)/V. Then we have, E(l) = V/(KA) We have 

complete the proof of identity (JU) 

By carefully analyzing, we find two properties of random walk model. One 
is the property of homogenous, that is, the properties, such as the locations and 
the associated probability, of the points are translation invariance. The other is 
the properties of isotropic, that is, the properties, such as the direction of next 
step and the associated probability, is rotation invariance. Therefore, we find 
the expectation of the length of random paths is the same as the expectation 
of the length of random chords. In order to demonstrate the fact, we use two 
dimensional model to show the result. In Figure 1, there is one random path, 
from point A to point B then from point B to point C. In Figure 2, there is 
one random path, from point D to point B then from point B to point E. The 
events of these two figures are equally likely. In the Figure 3, there are two 
random paths, one is the same as in Figure 1 and the other is the same as in 
Figure 2. In the Figure 4, there are two random paths, one is the path from 
point A to point B then from point B to point E. and the other is the path 
from point D to point B then from point B to point C. Using these facts, these 
two cases of events and their associated probabilities, events in Figure 3 and the 
events in Figure 4, must be equivalent to each other. These arguments support 
that the expectation of random path (Figure 3) is the same as that of random 
chords (Figure 4). 

The Theoretical Values of K 

Now, we are going to find the expectation of random chord. Without losing 
the generality, we solve the two dimensional problems instead of n-dimensional 
problem. Though there arc so called the Bertrand paradox in random chords 
problem, we are able solve paradox since the results computer simulation will 
show the correct answer of this paradox. This verifying work could not be done 
in that time. 

We consider a circle with radius r. The set of all chords is well defined. 
Therefore, we can find the expectation (average) of the lengths of the set of 
chords. Since the chords are uniformly and randomly distributed in the cir- 
cle, we can find some subset of these chords by defining a equivalent rela- 
tion. The parallel relation between chords is a equivalent relation on the set of 
that chords. The equivalent relation partitions this set into disjoint equivalent 
classes. Clearly, the expectations of these equivalent classes must be the same 
one. In the case of finite set, if the mean of all disjoint subsets is the same one, 
then the set has the same as that of these subsets. Therefore, we will compute 
the expectation of the random chords in one of these equivalent classes. In figure 
5, the measure of the distribution is dy/(2r). Therefore, we have expectation of 
random chords, E{C) 1 
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E(C) = f_ r 2xdy/(2r). (3) 



E(C) =irr 2 /(2r). 



(3) 



We interpretate equation ([3]) as E(C) — V n /V n -i since the diameter is the 
normal section of the circle. This interpretation can be generalize to higher 
dimension cases. 

From identity @, we have 



where V n -i is the volume of S"™ -1 the normal section of S n and A n is the bound- 
ary area of S n . Here, we mean that the sphere contains all its interior points. 
If it is necessary, we will use the term the sphere surface which is S n to be used 
in usual sense of mathematician. From equation ([J), we are able to calculate 
the theoretical values of all K. 

Discovery of some sequence 

The paper, by H. Chalkley, J. Conficld and H. Park, stresses how to estimate 
the ratio of volume and area. So do their followers'. But our interest is to find 
the values of K. It can be shown that the values of K depend on the dimension 
the space and that the values of K are not dependent on the shape or the size 
of the bounded regions. The ratio V/A of S n is well known. If we can find the 
expectation of random chords, then we are able to get the theoretical values of 



From equation (TTJ), We can estimate the value of by computer simulation 
while the theoretical value of K can be obtained by equation K = A/(VE(l)) . 
The values of are K x = 1/2, K 2 = l/ir, K 3 = 1/4, K 4 = 2/(3tt), ...K 7 = 5/32, 
...,^0 = 128/(31570, ... 

Our approach is to study the related topics rather than to solve this problem 
separately or independently. In order to convince the readers that the results 
we obtain are correct and can be verified by computer simulation, the computer 
programs have been run on PC, DEC VAX 9000, and IBM SP2 for more than 5 
years. It is computer that generates a phenomenon which we have never found. 
From this phenomenon, we try to find some mathematical model to fit it. It is 
quite natural to generalize the results to Riemannian manifold. And then the 
sequence 1/2, 1/tt, 1/4, 2/(3tt),. . . ,5/32,. . . , 128/(315tt),. . . , is discovered. So 
far, we have found that K^'s depend on the dimension of the space or manifold 
only. 

It takes time to judge how important the identity , E{1) — A/KV , is, be- 
cause there are many branches of mathematics such as geometry, analysis and 
probability or measure theory, concerning with this equation . Wc think that 
this is a interesting problem. 

We have done the computer simulation of random walk on Riemannian mani- 
fold, the surface of the sphere S n , n > 2. Almost the same result is obtained. 
We will organize another paper to show how to design the computer program 



K = K-i/A 



(4) 



K. 
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for simulating the random walk on the surface of sphere S n . We shall simply 
use Riemannian manifold S n in usual sense. In S 2 , we are shocked by a simple 
example. Let the boundary region be defined on the north pole of the earth 
S 2 . The region is a circle of which the center is the north pole and the radius 
R, R < nr/2, since R > nr/2, the region becomes a circle on the south pole. 
The chord of the circle is defined as the segment of a geodesic that contains at 
least one interior point of the circle, the bounded region. Of course, the two 
end points must be on the circle. Let us use the spherical coordinate system, 
p, 9 and <f>. Here p is constant since we are studying the problem on S 2 . For 
any 0, < < n/2, there is a circle defined on the north pole. The value of 
V is f® 2nr sin 6' 'rdd 1 , that is, 27rr 2 (l - cos0). The value of A is 2th- sin 0, The 
value of K is 1/n. Therefore, the value of E(l) or E{C) is nr(l — cos0)/ sin0, 
< < 7r/2. The circle is the north hemisphere of the earth when = ir/2. 
Then the value of Vis 27rr 2 , the value of A is 2nr and the value of E(l) or E{C) 
is nr. We find this is a circle of which all chords have the same length, nr. 
Since we were confused by this circle, we have spent more than one year's time 
to reinvestigate our formulation. Finally, we found the result is correct. 

References 

1. H. Solomon, Geometric Probability (Society for Industrial and Applied 
Mathematics, Philadelphia, Pennsylvania, 1978), pp. 1-32, 127-172 . 

2. H. Chalklcy, J. Cornfield, H. Park, A method for estimating volume- 
surface ratios, Science, 110, pp. 295-297 (1949). 

3. C. Kittel AND H. Kroemer, Thermal Physics (W.H. Freeman and Com- 
pany, San Francisco, ed. 2, 1980), pp. 399-414. 

4. R.P. Feynman et al., Lectures on Physics (California Institute of Tech- 
nology, 1963), vol. 1, section 43-3-section 43-5. 

5. K. Stowe , Introduction to Statistical Mechanics and Thermodynamics 
(John Wiley, New York, 1984), pp. 324-327. 

6. Hung-Ming Tsai et al., Estimation of the Mean Speed of Mosquitoes, 
Science Exhibition, Taichung, Taiwan (1993). 



5 



Figure 1: In this figure. There is a mosquito that gets into the region from 
boundary point A. Then the mosquito makes a turn at point B and goes out 
the region through the boundary point C. 
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Figure 2: In this figure. There is a mosquito that gets into the region from 
boundary point D. Then the mosquito makes a turn at point B and goes out 
the region through the boundary point E. 
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Figure 3: In this figure. There are two mosquitoes that get into the region, one 
of them gets into the region from boundary point A. Then the mosquito makes 
a turn at point B and goes out the region through the boundary point C. The 
other gets into the region from boundary point D. Then the mosquito makes a 
turn at point B and goes out the region through the boundary point E. 
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Figure 4: In the figure, there are two random chords, AE and DC in the circle. 
Clearly, The mean the lengths of these two chords is the same as that of the 
two random paths in Figure. 3. 
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Figure 5: In the figure, some random chords of a class are shown. The diameter 
that is perpendicular to the random chords is called the normal section of the 
circle. 
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the table of data 

dimension 1 dimension 2 dimension 3 dimension 4 dimension 5 
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The computer program 

program sphere 

! This main program, id and k are the control variable for 

! the dimension of the sace. m is the number 

! of partcles. n is the steps of random walk, 
dimension x(600100,10),dd(600100,10),dx(10),tow(10),tm(10),tn(10) 
open(2,filc=" D:\for\sp41. dat" ,status=" unknown") 
rewind(2) 

write(2,*)"t=0.01 pro=sp n=1000 4th aprali" 

pi=3. 1415926 

r=l 

t=0.01 

do id=l,10 

ii=id 

call brn(ii,a,ad) 

nn=a*1.15**id 

n=1000*nn 

m=1000*nn 

do k=l,10 

s=0 

ss=0 

! The m number of particles is uniformly distributed. 

! The partcles is generated uniforly in the general interval. 

! The method of rejetion is used to obtain the desired results, 
do i=l,m 

100 rs=0 
do j=l,id 

call random_number(xZ) 

x(i,j)=2*xz-l 

rs=x(i,j)**2+rs 

end do 

rl=sqrt(rs) 

if (rl .GE. r)go to 100 

end do 

do ll=l,n 

ssl=0 

! The directions of next step are generated, 
do i=l,m 
call rt(dx,id) 
do j=l,id 
dd(i,j)=dx(j) 
end do 
end do 

! The postions of m partcles are updated. 

! Of course, the counters s and ssl are used to 
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! count the number of partcles that hit the boundary. 
! The technique for using two counter s and ssl insure 
! the sum is correct. That is, sssl=ssl+l is not correc 
! counter, when the content is very large. 

do i=l,m 

do j=l,id 

dx(j)=dd(i,j) 

tow(j)=x(i,j) 

x(i,j)=x(i,j)+t*dx(j) 

tn(j)=x(i,j) 

end do 

rs=0 

do j=l,id 

rs=x(i,j)**2+rs 

end do 

rl=sqrt(rs) 

if(rl .LT. r) go to 300 

call tt(tow,tn,r,dx,id,t) 

do j=l,id 

x(i,j)=tn(j) 

dd(i,j)=dx(j) 

end do 

s=s+l 

ssl=ssl+l 
300 end do 

ss=ss+ssl 

ssl=0 

end do 

! The theoretical value of K, ty , is computed. 
ii=id 

call voa(ii,a,ad) 
tx=n 
sm=m 
b=a*id 

tk=ss/(sm*b*tx*t)*a 
ty=ad/b 

! The computer simulation of K is tk. 
write(2,*)tk,ty 
write(*,*)tk,ty 
end do 
end do 
close(2) 
stop 
end 

subroutine tt(tow,tn,r,dx,id,t) 

! In this subprogram, the new postion is determined 
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! when the partcle collides the boundary of sphere. 
! The elastic reflection is computed. A simple property is used. 
! That is, any position, including on the boundary, of the particle is vector 
which is 

! perpendicular to the (tangent of ) boundary 
dimension tow(10),tn(10),dx(10),f(10),g(10),h(10) 
ap=0 
pl=0 
do j=l,id 

ap=ap+dx(j)*tow(j) 
pl=pl+tow(j)**2 
end do 

tp=-ap+sqrt(ap**2+l-pl) 

pr=0 

prl=0 

do j=l,id 

f(j)=tow(j)+tp*dx(j) 

pr=pr+(tn(j)-f(j))*f(j) 

prl=prl+dx(j)*f(j) 

end do 

dl=0 

do j=l,id 

tn(j)=tn(j)-2*pr*f(j) 
dx(j)=dx(j)-2*prl*dx(j) 
dl=dl+dx(j)**2 
end do 

dl=1.0/sqrt(dl) 

do j=l,id 

dx(j)=dx(j)*dl 

end do 

return 

end 

subroutine rt(dx,id) 

! In this subprogram a random unit vector is generated for 

! direction of next step random walking, 
dimension dx(10) 

100 tl=0 
do j=l,id 

call random_number(xz) 
dx(j)=2*xz-l 
tl=tl+dx(j)**2 
end do 

if(tl .EQ. 0) go to 100 
tl=sqrt(tl) 

if(tl .GT. 1) go to 100 
tl=1.0/tl 
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do j=l,id 

dx(j)=dx(j)*tl 

end do 

return 

end 

subroutine voa(ii,a,ad) 

! In the subprogram, both the volumes of unit sphere and the normal section 
! section of the unit sphere are computed. For example, if input ii=3, 
! then output a is the volume of three dimensional unit sphere and 
! the other output ad is the area of a unit two dimensinal sphere, the disk. 

ii=ii-l 

call brn(ii,a,ad) 

ad=a 

ii=ii+l 

call brn(ii,a,ad) 

return 

end 

subroutine brn(ii,a,ad) 

! In this subprogram, the volume of unit sphere is computed 

! The parameter ii is the dimension, the parameter a is ouput and 

! the parameter ad is a redunduncy 

in=ii/2 

in=ii-2*in+l 

if (in .EQ. 2)go to 400 

f=l 

ls=ii/2 

do j= 1,1s 

f=f*j 

end do 

a=3.1415926**ls/f 
go to 500 

400 f=l 
do j=l,ii 
f=f*j 
end do 
ll=(ii-l)/2 
fl=l 
doj=l,ll 
fl=fl*j 
end do 

a=3.1415926**ll*2**ii*fl /f 
500 return 

end 
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